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Abstract 

We consider explicit and spontaneous CP-violation related to the profile of the two 
Higgs doublets at the MSSM electroweak phase transition. We find, in accordance with 
previous results, that in principle spontaneous CP-violation could exist in the MSSM at 
finite temperatures, but when the constraints from experiment and the strength of the 
transition are taken into account, the relevant region of the parameter space is rather 
restricted. Nevertheless, we show that in this small region, perturbative estimates 
need not be reliable, and non-perturbatively the region might be slightly larger (or 
smaller). To allow for more precise perturbative studies and for a lattice study of the 
non-perturbative corrections, we construct an effective 3d theory for the light stops 
and the two Higgses in the regime of large tuq. The 3d theory involves CP- violating 
parameters and allows to determine CP-odd observables. 
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1 Introduction 



The physics problem motivating the studies of the cosmological electroweak phase 
transition is that it could contribute to the matter-antimatter (baryon) asymmetry 
of the Universe [IJ. Several different ingredients are required for baryon asymmetry 
generation: one needs anomalous baryon number violating processes in the symmetric 
high temperature phase, microscopic C- and CP- violation, and thermal non-equilibrium 
(for a review, see 0). 

Recently, much progress has been made in understanding the order and strength of 
the electroweak phase transition (which determine whether there is non-equilibrium 
and what the rate of baryon number violation is in the broken phase), and also in 
understanding the rate of baryon number violation in the symmetric phase. It has 
been found that in the Standard Model, there is no electroweak phase transition at 
all for the experimentally allowed Higgs masses (for a review, see ||). In the MSSM, 
in contrast, there can be a first order phase transition, the strength of which has 
been studied perturbatively up to 2-loop level [Q-|§ and, because of the potentially 
bad convergence of the perturbative series, non-perturbatively employing dimensional 
reduction and three-dimensional (3d) lattice simulations P|— 1|12|| . The conclusion of 



these studies was that, from the point of view of the strength of the transition, Higgs 
masses up to 105. ..110 GeV are allowed for electroweak baryogenesis, provided that 
the right-handed stop is light, mi < m t0 p- 

As to the baryon number violation rate, its parametric form has been determined 
analytically (TT3], [TJ]], and numerical estimates exist both in the symmetric phase [IS] 



and in the broken phase [|16| of the theory. The latter case serves to determine the mag- 
nitude of higher order corrections to the perturbative 1-loop saddle point computation; 
the corrections were found to be as small as expected. 

However, it seems that the third requirement for baryon number generation, namely 
CP- violation, remains less well understood. If there were no CP-violation, one would 
produce the same amounts of baryons and anti-baryons, and no net asymmetry would 
arise. In the Standard Model, there is CP-violation as is experimentally observed in the 
/^-system, and this phenomenon is explained by the single complex parameter in the 
Kobayashi-Maskawa matrix. The Standard Model CP-violation is probably not large 
enough for the purpose of baryon number generation, though (for a review, see 0). 
Hence one would need new sources of CP-violation which, at the same time, do not 
produce any observable effects violating the existing experimental constraints. 

An interesting prospect for such a mechanism is provided by the MSSM. Indeed, 
there are new sources of CP-violation there. Fixing some regularization scale, one 
can redefine the fields so that the new CP-violating parameters appear in the trilinear 
couplings of the two Higgs doublets H%, H 2 and the squarks: 

£ C p = h t (A* t H 2 - (iHrfQJJa + H.c. , (1.1) 
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where H2 = iazH^, ht is the top Yukawa coupling, At, are complex parameters 
with the dimension of mass (A t is induced by soft supersymmetry breaking while /z is 
assumed to be the supersymmetric mass parameter appearing in the Higgs sector), Q a 
is the left-handed third generation squark, and U a is the right-handed third generation 
squark. However, if the CP-violating effects are proportional to the complex phases 
of At, fi, one is in the regime of explicit CP-violation, and then one again has to be 
careful in order not to violate the experimental constraints (for a review, see JF7| )[]. 

There is also another possibility for CP-violation in the MSSM. Indeed, since one 
has two Higgs doublets in the theory, there is the possibility of spontaneous CP- 
violation p|. This means that even though all the parameters were real, there could 
be a dynamically generated phase angle between Hi,H 2 (of course, without any ex- 
plicit CP- violation, an angle and minus the angle are equally likely). The possibility 
of spontaneous CP- violation is in principle there already at T = 0, but in practice the 
parameter values needed (in particular, a small CP-odd Higgs mass m^) are experi- 
mentally excluded [ 20|| . On the other hand, it may be that spontaneous CP- violation 



is more easily realized at finite temperatures [21 1, or even that it takes place only in 



the phase boundary between the symmetric and broken phases [22L I23I . A non-trivial 
CP-violating profile in the phase boundary could conceivably be quite useful for elec- 
troweak baryogenesis |24j-||31|, and at the same time, one would not need to worry 
about the constraints that have to be satisfied in the broken phase at T = 0. 



The parameter region found in [21, 22, ^3] for finite T spontaneous CP- violation 



consists of relatively small values of and relatively large values of tan/3 . Such values 
are not particularly ideal from the point of view of the strength of the transition, but 
could still be allowed in the small stop regime |J. 

In this paper, we analyze the general prospects for finite T spontaneous CP- violation 
in some detail. The emphasis is on understanding what kind of perturbative and non- 
perturbative effects there are, and whether the non-perturbative effects could be studied 
with lattice simulations. As a tool for these discussions, we construct an effective 3d 
theory which involves the CP-violating dynamics. Finally, we explain how this theory 
can be used for more precise perturbative and non-perturbative numerical studies than 
carried out here. 

It should be noted that all these finite T studies concern the thermodynamical equi- 
librium situation at T = T c . During the actual history of the electroweak phase transi- 
tion, the static equilibrium situation with T = T c is only reached if reheating to T c takes 
place after the transition, which is not generically true |55^ |S|. However, understanding 
the equilibrium situation is clearly a good starting point. 

The plan of the paper is the following. In Sec. we briefly describe the functional 
form of the relevant dimensionally reduced effective 3d theory, and the parametric mag- 
nitude of its couplings. In Sec. ^|we study, parametrically, what kinds of "instabilities" 

3 The experimental constraints might be weaker than often assumed, though |lq]. 
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could in principle occur. Such possible instabilities are the deviation of an effective 
tan/3 -parameter from a constant value within the phase boundary, spontaneous CP- 
violation, and the breaking of the charge neutrality of the vacuum. In Sec. |] we discuss 
the profile of tan/3 in some more detail, in particular with respect to whether it could be 
studied with lattice simulation. In Sec. [| we do the same for spontaneous CP- violation. 
We conclude in Sec. The details related to the derivation of the effective 3d theory 
described in Sec. || are presented in the Appendix. 



2 The effective 3d theory for H h H 2 ,U 

To discuss the physics of finite temperature CP-violation, we shall first construct the 
corresponding effective 3d theory, using the method of dimensional reduction applied 
to systems with phase transitions |33]| . The motivation for this approach is that it 
implements automatically the resummations needed at finite temperatures and is thus 
the simplest way of seeing which infrared (IR) problems can be cured and which remain 
non-perturbative. The 3d theory may also allow for lattice simulations, as discussed 
belowQ. 

We work in the light stop regime, so that the dynamical degrees of freedom appearing 
in the effective theory are the SU(2) and SU(3) gauge fields, the two Higgs doublets 
Hi, Hz, and the right-handed stop U (in principle, the U(l) gauge field should also 
be kept there, but it is expected to induce only small corrections |35[). Moreover, we 
assume that the left-handed stop mass parameter is relatively large, rriQ > 0.6 TeV, 
whereas the mixing parameters and the other mass parameters are not excessively so: 
\A t \ <rriQ, and and the other mass parameters are small compared with ttT. The 
details of the derivation of the effective theory are in the Appendix. We need here the 
main features only. 

The effective 3d Lagrangian is the most general gauge-invariant Lagrangian involving 
the given degrees of freedom: 



3d 



\p a p a + -G A G A 

4 ij v ^ 4 y v 



+ {Df H l )\D™ Hx) + (D?H 2 y(D?H 2 ) + {D1U) ] {D^U) 

+ mKTjHlH! + m 2 2 {T)HlH 2 + [m\ 2 (T)H{H 2 + H.c. ] + ml(T)lPU 

+ ryttfUHtHi + -i 2 U ] UH\H 2 + [ ll2 U^UHiH 2 + H.c. ] + \u{U^U) 2 

+ AiCF^tfi) 2 + \2{H ] 2 H 2 f + \iH\H x hIh 2 + XiHl^HlH-L 



+ 



\ 5 {H\H 2 f + X 6 HlH 1 HlH 2 + \ 7 H\H 2 H\H 2 + H.c. 



rt: 



t; 



rt, 



(2-1) 



It is interesting to note that the very good accuracy of dimensional reduction (in the Standard 
Model) has recently been confirmed also fully non-perturbatively with 4d finite temperature lattice 
simulations fl33l. 
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where Df = d { - igt a A a { , Df = $ - ig s T A Cf are the SU(2) and SU(3) covariant 
derivatives, t a = a a /2 where a a are the Pauli matrices, and U denotes the complex 
conjugate of the original right-handed stop field appearing in Eqs. (|1.1|) , ( |A.94j ). Of the 
parameters in Eq. (|2.1|), m 2 2 ,7 12 , A 5 , A 6 , A 7 can in principle be complex. The complex 
phases are due to the couplings A t ,\i in Eq. (|1.1| ). On the other hand, in the case of 
spontaneous CP-violation, all the parameters can be real. 

It should be noted that after dimensional reduction, the effective theory also con- 
tains the zero components of the gauge fields Aq, Cq, with mass terms (l/2)m\ AqAq, 
(l/2)rriQ CqCq, and, in principle, all the possible gauge-invariant quartic interactions 
involving Aq, Cq and the degrees of freedom in Eq. fl2.1|) (most of these can be con- 
structed out of the gauge-invariant bilinears AqAq, CqCq, H\Hi, H 2 H 2 , HiH 2 , WU). 
However, as argued in Appendix |A.5| , these terms are not essential for the present 
discussion, and thus Aq, Cq can be integrated out. 

It is also important to note that all other CP-violating operators, such as those 



related to FF (see, e.g., [24]), result in higher-order operators whose contributions 
to Higgs sector CP-violation are suppressed relative to the effects arising within the 
theory in Eq. Q J5§. 

The expressions for the couplings in Eq. ( |2.1| ) arising from dimensional reduction 
are summarized in Appendix |A. 7| (see also || |10|, [TT], |2l], |22], |2~3f). Here it is sufficient 
to know the leading parametric behaviour of the couplings: omitting multiplicative 
numerical factors and replacing logarithms by terms of order unity in the 1-loop terms 
proportional to l/(167r 2 ), we get 



^(T)~m?+(^ 2 -i^|A| 2 )T 2 , 
^(T)~m^+(^ + ^ 2 -^|i t | 2 )r 2 , 
m 2 12 (T) ~ m\ 2 + \h\k t ^ - J^Mtf + y^/A 
mUT) ~m 2 u+ (§<£ + \h\ - \h\{\Atf + |A| 2 ))t 2 , 



Ac/ ~ ^9s, 7i 
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~-hi\fi\\ l2 ~ht{l-\A 

9 2 hU m\ 2 M 2 /x* 
167r 2 [(2vrT) 2 (vrT) 2 



712 ~ h 2 t A* t fi* 



16tt 2 



A 3 ~^ 2 -</ 2 ), A 



1 2 



A, 



16tt 2 . 



4 m l2 

(2ttT)'- 



+ 9 



(M 2 /i 



*\2 
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A, 

A7 ^ 



1 

16vr 2 
1 

16tt 2 



-5 htAtfj, -g 



rn 



12 



(2ttT) : 



- 9 



-h$A* t j>* - g« 



ni 



12 



(2vrT) : 



- 



M 2/ i* 
(vrT) 2 
4 M 2 /i* 



(ttT)* 



(2.11) 
(2.12) 



where m 2 , m 2 , m 2 2 , m 2 / denote the zero-temperature mass parameters, M 2 is the SU(2) 
gaugino mass parameter, and 



A = 



A 

m Q ' 



(2.13) 



The signs (which do have some significance) have been obtained by assuming that itlq ~ 
1 TeV, T ~ 100 GeV, \A t \ ~ |//| < 100 GeV; for the general case, see Eqs. (|A.137|) - 
( |A.151 ). It should be noted that loop effects within the effective theory may generate 
contributions which are larger than those in Eqs. (|2.10] )- (|2.12|) ; see Sec. [|. 

We have written in Eqs. ( |2.2| )-( p7L2"D the couplings in 4d units, so that there is an 
overall factor Jq dr = 1/T in the action. This factor can, as usual, be defined away by 
a scaling of the fields, and the result is that the quartic couplings get simply multiplied 
by T. 

Note that the form of Eq. ( |2.1|) imposes several constraints on the parameters. There 
are various kinds of constraints (see, e.g., [0). First, the quartic part of the potential 
should grow (or remain the same) in every direction of the field space, for the field values 
that the effective theory is applicable. Otherwise the potential is "unbounded from 
below". Second, the standard electroweak minimum should be at least a metastable 
one at low temperatures. 

Denoting by t>i,t>2 the expectation values of the neutral components of Hi,H 2 , 
tan/3 = v 2 /vi, and 



7eff 



hi 



|/t| 2 cos 2 /? + (1 - I A I 2 ) sin 2 /? + (A* t (j* + A t fi) cos/3 sin/3 



(2.14) 



the first constraint says that for all values of j3 that lead to 7 e ff < 0, one must require 

1 



7eff< ^5^ 2 cos 2 2/3, 

where g 2 = g 2 + g' 2 . For instance, taking tan/3 = 
instability in the plane v 2 = unless 



h*\fi\ 4 < 



12 



(2.15) 

0, this means that there is an 

(2.16) 



Thus, we have to require at least that 0.4. There is a similar but weaker 

constraint for A t , A t < 1. 



5 



The second constraint applies in the broken electroweak phase, for the physical value 
of tan/3 . Then the simplest requirement is just that the right-handed stop mass squared 
be positive, which means that for the physical value of tan/3 , 7 eff must be positive (for 
the values of mfj we are interested in), and moreover, that 



2 

\A t - (x* cotp\ 2 > (2.17) 



m top 



This constraint is typically weaker than the one in Eq. (|2~T5|) . Another version of 
Eq. ( |2.17| ) applies at finite temperatures when one requires that the history of the 
Early Universe leads us to the standard electroweak minimum [|, 0, |). 
In order to satisfy these constraints, we assume that |A t | 2 , |/t| 2 <C 1. 



3 Parametric estimates for instabilities, 
and when is further reduction possible? 

Let us now inspect the Higgs doublet part of the action in Eq. ( |2.1| ). Taking into account 
the tree-level expressions for the parameters m^m 2 , (see Eqs. ( |A.137p , (|A.138|) ), we 
observe that the sum of the eigenvalues of the Higgs mass matrix is 

m\{T) + m 2 (T) = m 2 A + gs 2 + \h\ - \h\{\A t f + + 0^))t 2 ~ (gT) 2 . (3.1) 

Here vha is the CP-odd Higgs mass, we use the parametric convention g ~ h t ~ g,s, 
and we assume that \A t \ 2 + \fi\ 2 -C 1, see Sec. @. On the other hand, at the phase 
transition point, there is a direction along which one of the Higgs doublets is light, 



m 



2 

light 



{g 2 T) 2 (this is because the phase transition takes at the tree- level place when 
the second derivative of the effective potential vanishes in some direction). Then the 
other Higgs doublet must be heavy, m^ ea¥y ~ (gT) 2 , since the sum of the eigenvalues 
is constrained by Eq. (|3.1|). Thus the question arises, how could there be instabilities 
such as spontaneous CP-violation, since one of the Higgs doublets is heavy and could 
possibly be integrated out, as done in [|, |H], [TTJ] ? If there is just one dynamical Higgs 



doublet, there cannot be any spontaneous CP- violation. 

To establish rules for power-counting, recall that the first order electroweak phase 
transition is generated typically by gauge field (or stop) 1-loop terms of the form —Tm^ 
(or —Tm~). Moreover, let us assume that the quartic scalar self-coupling A is of the 
parametric order g 2 . Then, according to the standard argument, the phase transition 
takes place when the tree-level quadratic and quartic terms are of the same order of 
magnitude as the 1-loop contribution: 

4t-,2 2 3rri 3 2 4 /, r,\ 

g T v ~ g Tv ~gt>, (3.2) 
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where v denotes the light Higgs field expectation value. It follows that v ~ gT. Note 
that numerically, one would like to have in the broken phase that v>T due to the 
sphaleron erasure bound |], [|, but since the bosonic temperature scale means really 
~ 2itT, even these values can parametrically be thought of as ~ gT. The actual 
convergence of the perturbative expansion of course gets better with an increasing 
numerical value of v. 

Thus, we discuss values of order v ~ gT below. Let us stress that these estimates 
are assumed to hold also within the phase boundary, even though strictly speaking a 
fixed value of v is only obtained in one of the homogeneous (meta) stable phases. 



3.1 Estimates in the diagonalized theory 

Let us now diagonalize the two Higgs doublet model mass matrix (consisting of m\(T), 
m 2 {T), m\ 2 {T) in Eqs. ( p.2| )-( p^|) ) at the phase transition point. We denote m 2 = 
mf ight , M 2 = m^ eavy , and the corresponding fields by h, H, with vacuum expectation 
values Vh, vh (note that h, H are linear combinations of Hi, H2). Then, the diagonalized 
Higgs potential is of the form 

V{h,H) = m 2 h ] h + M 2 H^H + \i(h^h) 2 + \ 2 {H ] H) 2 + \ z h)hH ] H + \ i h ] HH ] h 
+ [\^H) 2 + \^hh ] H + \ 7 H ] Hh ] H + YL.c.]. (3.3) 

The couplings Ai, ...,A7 are some linear combinations of those in Eq. ( |2.l| ) (see, e.g., 
Eq. (6.21) in p[), but for simplicity we do not introduce a new notation. 

Since one would expect that the dynamics is dominated by the light field h, consider 
the set of Green's functions with only h in the outer legs. Conceivably one can, order 
by order in perturbation theory, organize the contributions of the heavy fields H so 
that they correspond to contributions to the parameters of an effective theory involving 
only h, 

C cS = (D™h)\D™h) + m^tf h + \ cS (tfh) 2 + ..., (3.4) 

where the terms not displayed are higher order operators. The expressions for the 
effective parameters m 2 s , X e g have been discussed in [f|: e.g., 

Aefi = Ai - ^(As + A3A4 + -Al + 2|A 5 | 2 + 12|A 6 | 2 - 6A*A 7 - 6A 6 A;). (3.5) 

The complication in the derivation of these parameters is that due to the interac- 
tions proportional to Xq, A7, which break the h — ► —h, H — > —H symmetries, there 
are reducible diagrams contributing to the effective parameters, related to the mixing 
between h,H generated at 1-loop level (see 0). 

However, clearly the construction of the theory in Eq. ( p.4[ ) cannot be valid for 
arbitrary quartic couplings: one must require that \{T jM -C 1 for the loop expansion 
parameters related to corrections such as those in Eq. (|3.5|) to be small. Parametrically, 
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this means that Aj < 0(g 2 ). Now, it is couplings at most of this order that arise from 



dimensional reduction, see Eqs. ( |2.6| )-( ]2TT2| ). Let us thus inspect whether this regime 



is consistent with instabilities such as spontaneous CP-violation. 

To analyse the kind of instabilities that can appear in the theory in Eq. Q2.1D , it 

is illuminating to write V(h, H) in terms of other sets of variables. Note first that 

V(h, H) can be written as a 2nd order polynomial in the explicitly real and gauge 
invariant operators 

O x = h)h, 2 = H ] H ) 3 = RetfH, A = \mh ] H. (3.6) 

These operators are constrained by the inequality 2 + 0\ < Oi<9 2 , but the existence of 
four operators nevertheless shows that there are, in general, four independent physical 
directions in the field space. However, the representation of the potential in terms of 
the Oj's is not exceedingly useful for studying these different directions, due to the 
constraint. 

For perturbative studies, a more useful representation is obtained by using global 
SU(2) and U(l) transformations to rotate h to unitary gauge and to remove one phase 
angle from H . Moreover, let us introduce an angular variable (3 by defining tan/3 = 
v 2 — v^+Vjj (note that this tan/3 is not the same as the usual zero-temperature 
parameter tan/3 = V2/V1 of the MSSM, since we have already made a rotation to a 
basis where m 2 2 (T) = 0, to arrive at Eq. (|3.3| )). Then, the two doublets can be written 
as 

v f cos/3 \ v ( sin (3 cos 9e l<t> \ 

k= V^\0 J' H= V^[sinPsin9 J" (3 ' 7) 



In terms of these variables, the tree-level potential in Eq. ( |3.3j ) becomes 

V(h,H) = ^ 2 (m 2 cos 2 /3 +M 2 sin 2 /3) 

+ ^v 4 { Ai cos 4 /3 + A 2 sin 4 /3 + ~ sin 2 2/3 [A3 + cos 2 6(\ A + 2A 5 cos2</>) 

+ sin2/3 cos6> cos</>(A 6 cos 2 /3 + A 7 sin 2 /3 )}. (3.8) 

There are now three possible types of deviations from the naive case that the heavy 
field H is zero and only the light field h is dynamical (i.e., j3 = 0). First, H can 
get a non-vanishing expectation value, corresponding to /3 7^ 0. As we will see, this 
deviation generically takes place and does not introduce any qualitative changes in the 
thermodynamical properties of the system. It turns out that (3 typically depends on 
the value of v. 

The other two changes are more drastic in the sense that they correspond to genuinely 
new "phases". Both of these deviations require that (3 7^ 0, see Eq. (|3.7|). The first 
possibility is that 9 can be non-zero. This means that the doublet H in Eq. ( |3.7|) has 
both upper and lower components non-vanishing. When gauge fields are taken into 
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account, this would mean that the photon becomes massive, and acts thus as a (non- 
local) order parameter for this phase (another way to say it is that the vacuum is not 
neutral). Clearly, the parameters are expected to be such that this possibility is not 
realized. 

The second possibility is that 0^0. This case corresponds to spontaneous CP- 
violation, if all the parameters appearing in Eq. ( |3.8| ) are real. Indeed, for the scalar 
fields, we can define the C-transformation as complex conjugation, while the P-trans- 
formation reverses the signs of the spatial coordinates but leaves the scalar fields intact 
(for a more general analysis, see f38|). Then 



h^h*, H^H*. (3.9) 



The action corresponding to Eq. ( [2.1|) is invariant under this transformation, for real 
parameters. However, the angle and the gauge-invariant operator 4 in Eq. ( |3.6|) are 
not invariant. Let us denote Ocp = O4. Ocp acts as an order parameter for this phase. 

The implications of (3 7^ 

Due to the fact that M 2 ~ (gT) 2 is large, we can assume all the three kinds of deviations 
to be small. For (3 this will be justified presently; for 9, <p this assumption is no 
restriction, since, without a loss of generality, we are choosing to inspect when the 
origin 6 = 0, <fi = becomes unstable (one may need to invert the signs of A6, A 7 for 
this). Expanding to second order in the angles in Eq. ( |3.8| ) and noting that Ajt> 4 ~ 
g 2 (gTf < M 2 v 2 ~ (gT) 4 , it is easy to see that 

(3 « sin/3 » ~ 0{g 2 ). (3.10) 

Here we assumed that A6 ~ C(5 ,2 ) ) f ~ M ~ gT. Thus, indeed, (3 is in general 
non-vanishing but small, and depends on v. 

Let us now inspect what a value (3 7^ means for the thermodynamics of the system 
and, in particular, for the construction of the effective theory in Eq. ( |3.4j ), which is 
supposed to account for the thermodynamics. As seen from Eq. Q3.8|) , a non-zero value 
of (3 causes a shift in the free energy of the order 

\m 2 v 2 ? 2 + \\ 6 v^ * ~ 0(g 8 T% (3.11) 



Is this contribution contained in the theory in Eq. (|3.4j) ? The answer is, no, since the 
term in Eq. ( |3.11| ) corresponds to a higher-order operator ~ (h^h) 3 . However, this 
contribution is suppressed by the relative amount 0(g 2 ) with respect to contributions 



within the effective theory in Eq. ( |3.4| ), since they are of the order 

m> 2 ~ A ef X ~ 0(g 6 T'). (3.12) 
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(The same estimate arises by comparing Eq. (|3.11|) directly with the tree-level 6-point 
function generated within the effective theory, ~ (^eff/ m eff) t ' 6 -) Due to the relative 
suppression 0(g 2 ), the negative sign of Eq. ( |3.11|) is also not a problem in the range 
of v considered. Thus, the effective theory is still useful and accurate in the phase 
transition region (for weak coupling), and (3 7^ does not change the properties of the 
system qualitatively. 

As a side remark, let us note that, for A6 ~ 0(g 2 ), the accuracy estimate 0(g 2 ) 
above differs from the one, 0(g 3 ), in |33], [36]. The reason is that with the present type 
of interactions, there can be reducible (left) diagrams with a heavy internal line, in 
additional to irreducible ones (right): 




In the Standard Model case considered in [^, in contrast, there are only irreducible 
diagrams, whose contribution is suppressed by one power of X 3 T/M ~ O(g) with 
respect to the reducible ones. In the MSSM, A6 is in fact typically parametrically 
smaller than 0(g 2 ), so that even there one may get errors smaller than 0(g 2 ). 



The implications of 0, 9 7^ 

Let us then look at the directions 0, 9. To quadratic order, the part of the effective 
potential depending on them is 

-v 4 {[3 2 [A 4 (l - 9 2 ) + 2Ab(1 - 9 2 - 2<p 2 )} + 2/3A 6 [l -\{9 2 + 2 )] }. (3.13) 



Thus, d> can be non-zero if 



/3 2 (4A 5 )+/?A 6 >0, 



(3.14) 



while 9 can be non-zero if 

/3 2 (A 4 + 2A 5 ) +f3\ 6 > 0. (3.15) 

Suppose that all the couplings Aj are of the same order of magnitude (~ g 2 ). Then, 
the (3 2 - terms in Eqs. (|3.14j) , ( |3.15|) are of order 0(g 6 ) since (3 ~ 0(g 2 ) and can be 
neglected, and the L.H.S. of Eq. ( |3.14| ) is, according to Eq. ( |3.10| ), 



/3A 6 



— — < 0. 

2 M 2 



(3.16) 



In other words, for couplings of this order of magnitude, the inequality in Eq. fl3.14j ) is 
never satisfied, is zero, and CP is not spontaneously violated! 
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When can then be non-zero? It is seen from Eq. Q3.14 ) that one must have A5 > 
and 

> « 5^ ~ < 3 - 17 > 

4p 2 ir 

(Recall that we assume v ~ M ~ gT.) Hence, CP can be spontaneously violated only 
if A5 is of the parametric order 0(1)! Clearly, in that case the construction of the 
effective theory in Eq. p.4j ) is not valid. 

Another possibility is that A7, which did not play any role in the previous estimates, 
is large. But to be effective, it has to be at least of the order 

A 7 ~ ^ ~ XeO(g- 4 ). (3.18) 

P 

If A6 ~ 0{g 2 ), A7 needs to be even larger than A5, and the reduction in Eq. ( |3.5|) again 
does not work. 

It is seen from Eq. ( |3.15| ) that the conditions for the breaking of 9 are similar as 



those for the breaking of 0, except that A4 < can have a further stabilizing effect. 

In conclusion, CP could in principle be spontaneously violated, and the photon could 
even become massive, if the quartic couplings of the two Higgs doublet model are 
suitable. In such a case, further reduction into a single Higgs doublet model does not 
work, because some of the expansion parameters 0{\T ' j 'M) would be large. However, 
this case is not a natural consequence of dimensional reduction in the MSSM, since the 
couplings produced this way are of too small a parametric magnitude, see Eqs. (|2lf)- 
( |2.12j ) (this statement is true even after radiative effects within the effective theory are 
taken into account, see Sec. |]). To get spontaneous CP-violation in practice, one thus 
needs somewhat large numerical factors and suitable parameter values, which overcome 
the parametric suppression. We will study these issues in Sec. [| 

3.2 Estimates in the original theory 

In the estimates so far we used for simplicity the diagonalized effective theory in 
Eq. (|3.3|). For completeness, it is perhaps good to see how the same estimates arise in 



terms of the original variables, before a redefinition which allowed us to get rid of m\ 2 - 
Let us start by reviewing the general tree-level condition for spontaneous CP-viola- 



tion within the theory in Eq. ( |2.1|) (this discussion corresponds to those in [21, P2j[). 
Putting now 6 — > 0, we choose, corresponding to Eq. ( |3.7| ), 

^ e -*)> H\H 2 = l - Vl v 2 e^. (3.19) 
Then, the part of the potential depending on <p is 

\^(cos 4>) = m\ 2 (T)v\V2 cos0 + - \§v\v\{2 cos 2 4> — 1) + X^v cos (ft + cos (ft 

(3.20) 
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Denoting again tan/? = 1)2/ V\, we note that the second order polynomial V^cos^) has 
a non-trivial minimum, if A5 > and 



l/GMI 



m\ 2 {T) + (1/2) (A 6 cos 2 /3 + A 7 sin 2 /r 



2A5 sin/3 cos/3 u 2 



< 1, 



(3.21) 



-f((3,v). These are the standard (tree-level) conditions for the 



in which case cos 1 

existence of spontaneous CP-violation [|1^ . 

Suppose that A6, A7 < 0(g 2 ), as one would like to have. Then, 
order 0(g 4 T 2 ) to compensate for the terms proportional to X e v 2 , \ 7 v 2 in Eq. (|3.21|) . 



nr 



12 



(T) must be of the 



Since m 2 2 (T) is small, one of the masses m 2 (T),77?|(T) has to be large, ~ {9T) 2 , in 



order for the eigenvalues to sum up to Eq. (|3.1| ). Let us assume^] that this role is played 
by : m 2(^) = M 2 + 0(g 2 T) 2 . Then, it can be seen from the Higgs part in Eq. ( |2.1| ) 
that 

ih ^_<+^iv soW (322) 



But this means that the term proportional to A7 in Eq. (|3.21|) can again be neglected, 
unless A7 is very large, and that Eq. ( |3.21|) becomes 



\f(M\ 



M 2 sin/3 



M 2 



2A 5 sin/3 v 2 2X 5 v 2 ' 



(3.23) 



Thus, we are again lead to the parametric estimate A5 > 0(1) in order to have 
\f(/3,v)\ < 1. Large numerical values of v (a strong transition) allow for smaller 
numerical values of A5. 

To summarize, we have found that according to tree-level considerations, it is not 
enough to have A 5 > to get spontaneous CP- violation, but one should have A 5 > Oil). 
The reason is that with the type of mass parameters and couplings that appear in the 
MSSM, the numerator and the denominator in Eq. (|3.21| ) tend to vanish at the same 
point, when the equations of motion for sin/3 are taken into account. 



4 The wall profile 

In the previous section, we saw that even the heavy Higgs doublet H can have a 
non-zero (even though small) expectation value in the broken phase, and a corre- 
sponding non-trivial profile at the phase boundary (since (3 typically depends on v, 
see Eq. (|3.10|) ). In the usual perturbative language in terms of the doublets Hi,H 2 , 
the possible non-trivial profile means that the ratio of field values, tan/? = 1)2 /v%, 
can change within the phase boundary. Since the typical values of H are very small, 

5 For completeness, let us note that in the MSSM it is actually m^(T) which is large, but here it is 



convenient to exchange the fields, to keep the discussion as close to Sec. 3.1 as possible 
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(H) ~ (h) tan/3 ~ g 3 T (see Eqs. (|3.2j ), (|3.10|) ), tan/5 is expected to change only by a 
small amount. Nevertheless, in P7|, |29[ it is argued that the deviations of tan/3 from a 
constant value are important for baryon number generation in the MSSM (this may not 
be a necessary requirement, though The purpose of this section is to discuss 

briefly the prospects for perturbative and non-perturbative determinations of tan/3 . 

Note first that, parametrically, the profile of tan/? , or the expectation value of H, is 
not computable in perturbation theory. This is simply because, as can be seen from 
Eq. (|3.10|) , the value of H is determined by that of h, which in turn is non-perturbative 
for h ~ gT (i.e., at least within the phase boundary). Numerically, the convergence 
might of course happen to be reasonably good, especially if the transition is strong as 
is required for baryogenesis. The wall profile obtained from the 2-loop Landau-gauge 
effective potential by solving the classical equations of motion has been computed 
in H (see also g, gjj). 

Let us then discuss whether the profile of tan/3 could be determined with lattice 
simulations. We first give a gauge-invariant generalization for what one may mean 
with the profile of tan/? . A static profile can only exist at the critical temperature. 
Then one may define 

\ v 1 = (Hj Hi) — (H-Hi) , tan 2 /3 = 4, (4-1) 

2 symm. yf 

where i — 1,2 and the latter expectation value is taken in the homogeneous "symmet- 
ric" high-temperature phase. This definition is clearly gauge and scale independent. 

Suppose now that we have a lattice geometry such that planar interfaces are favoured. 
There is a zero mode related to the location of a planar interface, but this can in 
principle be removed. Averaging over the set of configurations and over the (£1,2:2)- 
plane, we get a value for tan/3 (£3). This definition should not be sensitive to the lattice 
spacing. However, this definition is sensitive to the cross-sectional area of the lattice. 

Indeed, denote by A(x), x = (xi,x 2 ), the deviation of the interface location from 
the average. Then a configuration with A(x) 7^ has the action 

S- mt ^a J d 2 x^[d t A{*)]\ (4.2) 

where a is the surface tension. The fluctuations of the measured projection of the 
interface location A are distributed as p(A) = / d 2 x(5(A(x) — A)). As a consequence, 
the effective width I of the interface seen after averaging over the area A, is expected 
to diverge as 

2 J dAA 2 p(A) 1 r 2 2 r 2, 1 1, A 

where £ is of the order of the longest Higgs correlation length. 
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In conclusion, the profile of the interface can only be studied by averaging over a 
finite area. This can be accomplished either by restricting the total cross-sectional 
area, or by averaging the surface profile only over a finite subarea of the large total 
area ('coarse-graining'). However, we do not regard this as a fundamental problem. 
From the practical point of view, a more serious concern is that since the deviations 
from a constant value of tan/3 are expected to be very small, they are most probably 
not visible from below the statistical noise. 

5 Spontaneous CP-violation 

In Sec. |3]we saw that, parametrically, spontaneous CP- violation is not naturally realized 
in the MSSM. However, parametric estimates need not always be reliable due to the 
actual numerical factors. In this section, we thus look at the issue from a more practical 
point of view. We will first make some perturbative estimates and then comment on 
the possibility of lattice simulations. 

Let us start by formulating the problem as follows. The question concerning the 
existence of spontaneous CP-violation can be factorized into two parts: 

(1) The perturbative computation of the parameters of the action in Eq. fl2.1|) . This 
computation is not sensitive to the infrared problems of finite temperature field 
theory and is thus in principle well convergent. 

(2) The determination of the constraint that the parameters of the 3d theory have to 
satisfy in order for there to be spontaneous CP- violation. As we will discuss, this 
constraint is sensitive to the infrared problems of finite temperature field theory 
and is ultimately to be determined non-perturbatively. As a concrete example 
of the non-perturbative nature of the problem, note that the parameter m 2 2 (T) 
has a logarithmic scale dependence on the 3d renormalization scale parameter A, 
m i2(^) ~ (71712 + ...)T 2 /(167r 2 ) ln[A/(g 2 T)], and the relevant scale for, say, the 
statement m 2 2 (T) = 0, cannot be fixed perturbatively. 

Concerning item (1), we will assume the expressions in Eqs. ( |A.137|) -( |A.151| ) to be 
precise enough for our purposes. If needed, the accuracy can be improved upon on 
two aspects, in particular: the zero-temperature MS-scheme parameter m 2 2 (/i) could 
be fixed even more precisely in terms of the physical observables of the theory by a 
more accurate vacuum renormalization, and the finite temperature corrections to the 
parameter m 2 2 (T) of the effective 3d theory could be derived up to 2-loop level, to fix 
its scale dependence (see, e.g., [^3| and the last Ref. in fllTf). 

Concerning item (2), we will here work basically at tree- level, estimating the effects 
of loop corrections only in the regime where they can be reduced to the construction 
of a further simplified effective theory. We argue that this approximation contains the 
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largest uncertainties in the present computation, and we explain in Sec. |6] how it can 
be improved upon. 



5.1 Tree-level constraints inside the 3d theory 



The (first) tree- level condition for spontaneous CP- violation is given in Eq. ( |3.21| ). To 
see when it could be satisfied in practice, note first that for very large fields v , 



l/GMI 



A 6 cos 2 (3 + A7 sin 2 /? 



2A 5 sin2/? 



>1, 



(5.1) 



since typically |A 5 | <C |A 6 |, |A 7 |, see Eqs. ( [2.10|) -( p.l2|) . For very small fields v — > 0, on 
the other hand, 



\mv)\ 



m 



12 



> 1. 



(5.2) 



A 5 sin2(3 v 2 

Moreover, the function (a + bv 2 )/ (cv 2 ) is a monotonous function of v 2 . Thus there can 
only be a region with \f(/3,v)\ < 1 provided that the numerator of Eq. ( |3.21| ) is zero 
at some field values. This requires that 



|m? 2 (T)|~{|A 6 |,|A 7 |K 

and that the signQ of m 2 2 (T) is opposite to that of A6 cos 2 /3 + A7 sin 2 /3 . 
The second tree- level constraint concerns A5: 

A B > 0. 



(5.3) 



(5.4) 



The larger A5 is, the larger is the range of values of v where spontaneous CP-violation 
can take place. In Sec. ^T|] we got a stronger constraint, A 5 > 0(g 2 ) (assuming now 
that v ~ T in Eq. ( |3.23| ), instead of v ~ gT as in Sec. |3.2| ). However, since we will 
see that the uncertainties in the constraint related to A5 are large (larger than those in 
the constraint related to m 2 2 (T) in Eq. ( p.3| )), we consider here only the most relaxed 
version of the As-constraint, the one given in Eq. (|5.4j ). The true requirement for A5 is 
eventually to be determined numerically. 

Note now that, according to the expressions for the parameters related to CP- 
violation (m 2 2 , A 5 , A 6 , A 7 in Eqs. (|2~1 ), (|2.10|) - (|2.12|) ), the inequality in Eq. ( |5~^) can be 
satisfied. However, assuming that v< T, we see from Eqs. ( |2.11| ), ( |2.12| ) that Eq. ( |5.3|) 
goes over into 
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(T)\< 



[J4_ 

1167T 1 



\A* t /l*\T 2 , 



g* \M 2l x* 



16tt 2 



7T- 



(5.5) 



6 Since we are in this section considering spontaneous CP-violation, we assume that At, jl are real, 
but we nevertheless keep the same notation for them which arises in the complex case. 
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If the expression for m 2 2 (T) in Eq. (|2.4j) is dominated by contributions from some single 
particle species, then Eq. ( |5.5| ) clearly cannot be satisfied, as all the different types of 
terms in Eq. ( |2.4|) are much larger than needed for Eq. ( |5.5|) . This holds independent of 
whether the dominant effects come from squarks (A t , (x 7^ 0), gauginos and Higgsinos 
(M 2 , n ^ 0), or Higgses (A t = // = 0). 

Thus, the only way to get spontaneous CP- violation is to have a cancellation between 
the contributions of the different particle species. Indeed, one can get m 2 2 (T) ~ at 
relatively large values of rriA, if sign(Aj/i*) > 0: according to Eqs. (|2.4|) , (|A.139| ). 

m 2 12 (T) ~ -\m\ sin2/3 + \h 2 t A* t (l*T 2 . (5.6) 

This behaviour is obtained for mfj <C (2ttT) 2 <C tuq. Numerically, Eq. (|5.6j ) could be 
close to zero if, for instance, T ~ T c ~ 85 GeV, tan/3 ~ 15, A t ~ 0.6, ft ~ 0.3, and 
77u ~ 70 GeV.Q 



5.2 Stop loops inside the 3d theory 

In the discussion so far, we have completely ignored loop effects within the effective 
theory in Eq. (|2.1| ), which might change the estimates noticeably (indeed, it is only loop 
effects which induce a first order transition and a stable phase boundary in the first 
place). We now proceed to estimate the loop effects related to the stop field U in the 
case that it is heavy enough to allow for a perturbative treatment. By computing its 
effects on the parameters of the Higgs sector and applying then the tree-level estimates 
in Eqs. ( |5.3| ), ( |5.4| ), we show that the loop effects can be large in the regime where 
m 2 2 (T) given in Eq. (|5.6|) is small. On the other hand, if one is clearly outside this 
regime, then loop effects cannot trigger spontaneous CP-violation. 

To be in the perturbative regime, we assume that the stop mass mfj(T) is of the 
parametric order (gsT) 2 . In other words, these estimates do not apply when one is 
very close to the "triple" point where the direction U can get broken (see |7], §]). 

Inside the 3d theory, 1-loop contributions to the masses are of the form 

d 3 p 1 m 

(5.7) 



(27r)V + m 2 4vr' 
while 1-loop contributions to quartic couplings are of the form 

d 3 p 1 1 



(27r) 3 [p 2 + m 2 ] [p 2 + m 2 ] An {m a + m h ) 



(5.8) 



7 It is interesting to note that the negative tree- level term in m\ 2 (T) can also be compensated for 
by a large stop vev ~ h^A^fi* (WU) in the stop breaking phase ||, However, there the Higgses 
are in the symmetric phase. 
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The stop contributions are 



5m\ 2 (T) = -^- ll2 Tmu{T) ~ -h 2 A* t fi*g s T 2 } (5.9) 

6X5 = 777a ~ 777 , (5-10) 

167T mu{T) g s T 

%nmu(T) mmu(T) 

We observe that the contribution to the mass parameter m\ 2 (T) is of order 0(g 3 )T 2 
and is parametrically small compared with the leading T 2 -term. The contributions 
to the quartic couplings are, in contrast, of order 0(g 3 ) and thus larger than the 
corresponding tree-level terms in Eqs. (|2.10| )- (|2.12| ). However, even these corrections 
are not large enough to make Eq. ( p.3| ) naturally satisfied for v ~ T, unless the zero 
temperature and g 2 T 2 -terms in m 2 2 (T) cancel, which was the requirement in Eq. (|5.6|). 

The negative stop contribution to A5, in fact, tends to make spontaneous CP- 
violation less likely, since Eq. ( |5.4j) is more difficult to satisfy. This effect should be 
compensated for by the gaugino-Higgsino term in Eq. ( 2.10|) . 



Finally, consider what happens in the limit that the stop contributions become non- 
perturbative, mjj(T) ~ g\{I '. In this case, one cannot look at contributions to effective 
parameters, but one may try to estimate the effects numerically by summing the con- 
tributions to the effective potential. For real parameters (A* t — > A t , fi* — > jl), the 
1-loop stop contribution to the effective potential is of the form 



V ~ — — 

2vr 



1 1 - I 3 / 2 

m 2 v {T) - -h 2 t v 2 Afi\ 2 + -h 2 v 2 2 (l - \A t \ 2 ) + h 2 Vl v 2 A t fx cos ^ . (5.12) 



For large mfjiT) and sign(A£/t*) > 0, this term is minimized by cos0 = 1, which 
again corresponds to no CP- violation. However, for small m^T) the term becomes 
non-analytic and may induce non-trivial behaviour. At the same time, higher order 
corrections are of the same magnitude as the 1-loop term. 

To summarize the effects of the stop loops: m 2 2 (T) cannot be made small unless it 
already is so, and thus spontaneous CP-violation cannot be triggered far away from 
where it was expected to take place according to Eq. ( |5.6|) . Moreover, even in this 
regime, stop loops in fact seem to make spontaneous CP-violation less likely at least 
as long as they allow for a perturbative treatment, since they give a large negative 
contribution to A 5 . 

5.3 Higgs loops inside the 3d theory 

Although the Higgs loops involve effects from the non-perturbative mass scale ~ g 2 T, 
it is nevertheless illuminating to inspect what kind of contributions they would give to 
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the effective parameters according to the loop expansion. Note that the gauge fields 
(and Aq,Cq), in contrast, do not contribute to the CP- violating couplings at leading 
order. 

It turns out that the Higgs field contributions to the CP-violating parameters are 
typically smaller than the stop contributions. For instance, in the limit m 2 2 (T) « 0, 
there must be at least one appearance of the couplings A5...A7, which are of order 
hf/ (16tt 2 ) instead of order h 2 as the couplings 71,72,712: 



5X5 _ (5/2)A| + 2AiA 5 (5/2)Aj + 2A 2 A 5 6A 4 A 5 + 4A 3 A 5 + A 6 A 7 
T 87rmi(T) 87rm 2 (T) 4vr[mi(T) + m 2 (T)] ' 

<5A 6 12AiA 6 3A 3 A 7 + 2A 4 A 7 + 2A 5 A 7 3A 3 A 6 + 4A 4 A 6 + 10A 5 A; 
~T" _ 87rmi(T) 8nm 2 (T) 4vr[mi(T) + m 2 (T)] ; 

5X 7 _3A 3 A 6 + 2A 4 A 6 + 2A 5 A* ^ 12A 2 A 7 ^ 3A 3 A 7 + 4A 4 A 7 + 10A 5 A 7 
~T~ ~ 87rmi(T) 87rm 2 (T) 47r[m 1 (T) + m 2 (T)] 



(5.13) 
(5.14) 
(5.15) 



The terms where there appears only the light Higgs with mass ~ g 2 T are seen to 
contain non-perturbative contributions which are parametrically as large as the original 
couplings. There are also other types of contributions of similar magnitude: expanding 
formally in m\ 2 (T) 7^ 0, one gets for example terms of the type 

^ *mm (5 is) 

where A ~ {Ai,...,A 4 }, i = 1,2. All these contributions involve expansion parameters 
~ {Ai, A 4 }T/((yf 2 T) ~ 1, and are thus not perturbatively computable. 

5.4 Numerical values 

Let us finally look at the numerical requirements for the parameter values after taking 
into account the dimensionally reduced effective couplings in Sec. |A.7| , together with 



the stop contributions in the regime where they are still assumed to be perturbative, 
Sec. [5T2"| . To compensate for the negative stop contribution to A5 in Eq. (|5.10|) by a 
positive gaugino-Higgsino contribution in Eq. ( |A.149|) , one needs that 



• rriQ is large, because then the gaugino-Higgsino term suppressed only by temper- 
ature is relatively more significant; 

• M 2 is large in order to increase the gaugino-Higgsino contribution, but A t is not 
very large, to decrease the stop contribution. 

On the other hand, to allow for relatively large values of rriA at the point where m 2 2 (T) 
in Eq. (|A.139 ) is small (for fixed T ~ T c ), it is preferable to have sign(M 2y u*) < and 
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large (this may help also in making A6,A7 positive, so that m\ 2 (T) may remain 
slightly negative). 

Using the expressions in Sec. |A.7| and those in Eqs. ( |5.9|) , ( p-lOp , we estimate that 



numerically the best region for these requirements is 

m Q >800 GeV, M 2 > 70 GeV, |/i|>50GeV, |A|<0.2, (5.17) 

where A t is as defined in Eq. (PTT3). In these estimates, we have assumed that T ~ 
T c ~ 85 GeV, as can be observed by equating Eq. ( |A.138|) with zero in the limit that 
tan/3 is large (for numerical estimates of T c based on the 2-loop potential and lattice 
see, e.g., |6|, [D], fH|). Nevertheless, even in the regime of Eq. ( |5.17| ) there remains the 
constraint 

m\ sin2/3 < 400 GeV 2 . (5.18) 

Hence one needs, say, itla ^65 GeV and tan/5 > 20. Clearly this parameter region is 
at most barely consistent with what is required for the strength of the transition in 
order to preserve the baryon asymmetry generated, and with experimental constraints, 
which both provide a lower bound for ||. 

5.5 Prospects for lattice simulations 

In the case of explicit CP- violation, some of the parameters appearing in the action have 
complex phases and the expectation values of CP-violating operators such as O4 = Oqp 
in Eq. ( |3.6| ) are proportional to these. Due to the experimental constraints on the 
complex phases which are thought to be relatively strong, the magnitude of explicit 
CP-violation should be relatively small for realistic parameter values. In this case, a 
lattice study of the CP- violating effects appears rather difficult from the practical point 
of view, for the same reasons as discussed in Sec. ^ for tan/5 . Perhaps the best that 
could be done is to solve the classical equations of motion for the CP-violating wall 
profile using the 2-loop effective potential. Even though the effects may be small, this 
case is of interest since it has been argued that explicit CP-violation might suffice for 
baryogenesis [^9], Note also that large explicit effects may not be excluded in 
which case also a lattice study is feasible. 

In the case of spontaneous CP-violation, on the other hand, the parameters can be 
real and yet the effects can be large and well detectable: the distribution of Oqp should 
have two peaks at opposite values of Oqp (perhaps only within the phase boundary), 
which gives an unambiguous signal. The role of explicit CP-violation is merely to 
choose one of the peaks as the physically preferred one. 
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6 Conclusions 



In this paper, we have studied the general prospects for spontaneous CP-violation 
around the electroweak phase transition in the MSSM. Parametrically, spontaneous 
CP- violation seems not to be easily realized. Nevertheless, we have shown that numer- 
ically there is a small regime, given in Eqs. ( fS.lTQ , Q5.18 ), which could be favourable for 
spontaneous CP-violation. We have shown that in this regime, the IR-sensitive loop 
contributions related to the stops are large and even non-perturbative in the light stop 
regime. The Higgs field contributions are non-perturbative, as well, but are numeri- 
cally smaller than the stop contributions. However, it does not seem obvious that the 
large effects would naturally increase the likelihood of spontaneous CP-violation. It 
is rather that one needs delicate cooperation from different particle species to get an 
effect. Moreover, there is the concern that the regime found for tua within the present 
approximations, is almost certainly excluded due to experimental constraints and the 
requirement for a strong 1st order phase transition ||. 

It is conceivable that in many extensions of the MSSM, such as the NMSSM and 
various two Higgs doublet models, the parameter region allowing for spontaneous CP- 
violation would be larger than in the MSSM. Moreover, even if there is no spontaneous 
CP- violation, it could be that in the case of explicit CP-violation, the non-perturbative 
dynamics of the 3d theory affects strongly the explicit effects (especially if the explicit 
effects are large, as proposed in |18|). It has also been argued that there could be CP- 
violating solitons in the broken phase of the theory, even though there is no spontaneous 
CP- violation [4~2~|. Thus, in spite of the fact that the parameter space for spontaneous 
CP-violation in the MSSM seems to be quite small, we believe that a more precise 
study of these questions is motivated. We expect relatively large uncertainties in the 
present results. 

The present results were based on deriving an effective action for the IR degrees 
of freedom of the theory, and studying then this effective action at the tree-level. 
The first part thereof, the derivation of the effective action, is a purely perturbative 
computation, but even its accuracy can in principle be improved upon. In particular, 
the MS parameter ml 2 (p,) could be related more precisely to physical observables such 
as niA, and the parameter m\ 2 (T) of the effective 3d theory could be derived with 
2-loop accuracy. Nevertheless, we expect that it is the non-perturbative IR-properties 
of the theory that are responsible for the largest uncertainties. 

The first step of a more precise IR-study would be to compute the full 1-loop and 
2-loop effective potentials within the effective theory, and then to solve the complete 
classical equations of motion for the wall profile consisting of v\,v 2 , <p, to see when <p can 
deviate from zero. Solving the classical equations of motion is in principle straightfor- 
ward but can lead, in practice, to somewhat difficult instabilities and slow convergence; 



a working algorithm has been recently proposed in [40 



At the second step, the theory in Eq. ( |2.1| ) could be studied with lattice simulations. 
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Concerning the lattice study, it should be realized that the wall profile can only be 
defined with a finite cross-sectional area, see Sec. |. However, we do not consider this 
to be a fundamental problem (the area dependence is only logarithmic). It should be 
noted that a lattice study of spontaneous CP-violation is, from a practical point of 
view, more promising than for instance a study of the non-trivial profile of tan/3 , since 
the effect can be observed by finding a non-trivial (two-peak) distribution of a single 
operator, such as Ocp in Eq. (|3.6|). The breaking of global discrete symmetries in 3d 



gauge theories has previously been studied numerically in for the SU(3)+adjoint 
Higgs model. 

Finally, it is interesting to note that the leading parity violating effects are very dif- 
ficult to study with lattice simulations, since they are induced by higher-dimensional 
operators in the 3d theory which are complex in Euclidian space ]36| . The C- violating 
(and thus CP-violating) effects discussed in this paper are, in contrast, contained in a 
real super-renormalizable Euclidian action, and there should be no fundamental prob- 
lems in the simulations. 
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Appendix A 

In this Appendix, we derive the 1-loop expressions for the parameters in Eq. ( [2.1Q . 
CP-violation is assumed to appear in the original theory (at some renormalization 
scale Jxq) only in the trilinear couplings A t ,fi, defined in terms of the Euclidian 4d 
Lagrangian as in Eq. (|1 . 1|) . Thus, in particular, the Higgs sector mass parameter 
m\ 2 and the gaugino mass parameter M 2 are assumed to have been tuned real at jl , 
employing field redefinitions. Of the parameters in Eq. ( |1.1| ), fi is taken to be equal to 
the supersymmetric mass parameter which appears also in the Higgsino mass matrix, 
although in principle there could be an arbitrary soft component in the trilinear term 
to which /i contributes pH| . The soft parameter A t , in contrast, only appears in the 
interaction in Eq. ( |1 . 1| ) . 

We shall work in the limit m 2 <C (2itT) 2 <C tuq, A\ -C mg, where m 2 signals all the 
mass parameters of the theory other than m 2 Q. m 2 = {m 2 , m 2 , m 2 2 , mf/, |/^| 2 , M 2 , M 2 }. 
Note that the trilinear coupling A t need not, in principle, be smaller than 27rT. 

In this limit, the construction of the effective theory can be carried out in two steps. 
First, in Sec. |A.1|, we integrate out the left-handed squark field Q at zero temperature. 



Second, we go to finite temperature and integrate out the non-zero Matsubara modes 
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of the remaining bosons (Sec. |A.2|) , as well as all the fermions (Sec. |A.3|) . For the sake 



of brevity, we discuss the field redefinitions related to all these steps at once in Sec. [A. 4 



In Sec. [A.5| we discuss issues which are somewhat less central to the present problem, 
namely the values of the gauge couplings and the removal of the zero Matsubara modes 
of the gauge fields A ,C . Finally, for completeness, we discuss the most important 
effects related to vacuum renormalization in Sec. IA.6I. All the results are collected 



together and simplified final expressions for the parameters appearing in Eq. ( |2.1| ) are 
given in Sec. |A.7| . 

It should be mentioned that many of the expressions given here have also been derived 
in other places. In particular, the dimensional reduction steps (but for tuq <C (2ttT) 2 
and mostly only for the Higgs sector) have been carried out in [P[-[|IT] (see also [|45|j ). 
In addition, expressions for the quadratic and quartic Higgs couplings involving CP- 



violation have been derived, in various limits, in [21, f2j[. 

The notation for the parameters of the final effective action, to be used below, is 
given in Eq. fl2.1p . In addition to the effective parameters, the individual reduction 
steps also give contributions to the wave functions, which we denote by 

5C = {D™H 1 )\D™H 1 )8Z 1 + {D™H 2 )\D™H 2 )5Z 2 + (£>?£/) t (£>•£/) ^ 

+ [{Df^yiDfH^Zu + R.c.]. (A.l) 

We do not write down the MS scheme divergences 1/e in any of the expressions; these 
can be restored by In ft 2 — > 1/e + In ft 2 . In the loop corrections, we always neglect terms 
proportional to the U(l) gauge coupling g' 2 . 

A.l Integrating out Q at T = 

Let us consider the case that the mass parameter related to the left-handed stops is 
much larger than the temperature, uiq 3> 2nT. Then, the field Q can be integrated out 
from the zero temperature theory. Note that this case is different from the consideration 
in [46|, where all the squarks were integrated out with degenerate mass parameters. 
Evaluating some typical finite temperature integrals numerically, we expect the results 
to be well reliable when, say, ttiq > 0.6 TeV (T c < 100 GeV). 

The left-handed stop field Q contributes to the effective action of the two Higgses 
Hi, if 2 and the right-handed stop U already at tree- level, due to graphs of the type 

s O ' N ' 

* s / \ 

/ N ✓ \ 

where the dashed line denotes the light degrees of freedom, the solid line the heavy 
field Q, and the filled circle an effective vertex. 
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Using the notation A t , /t denned in Eq. (p,13|) , the tree-level couplings of the action 
in Eq. (p.l|), after integrating out Q, are 



9l 



Ac/ = ^, 7i = -/i t V, 7 2 = ^(1-^), 7l2 = /i 2 A*/r 
o 

> > _ g 2 + £ /2 > g 2 - g' 2 , g 2 

Ai — A 2 - , A 3 - , A 4 - -y, 

A5 — A6 = A7 = 0, 



(A.2) 



where /it is the top Yukawa coupling, and gs, g, g' are the strong, weak and U(l) gauge 
couplings. There are no contributions to the wave functions in Eq. (|A.1| ) at tree-level. 

At 1-loop level, the following types of contributions arise (we expand in m 2 /mQ, 
q 2 /m 2 Q ): 



Q. 
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fflg 167T 2 ' 



(A.7) 
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1 



(p 2 + rriQ) 2 (p 2 + m 2 ) 2 
1 1 



dp 



(p 2 + m 2 ) 2 
(A.8) 



where dp = <i 4_2e p/(27r) 4_2<E . Here the subtracted parts are contributions within the 
effective theory: after the subtractions the results are IR-safe as they should be. The 
effective six-point vertex in Eq. ( |A.7| ) is shown just for illustration; it gives a contribu- 
tion suppressed by the relative amount m 2 /m 2 Q and is not kept in the actual theory in 
Eq. Q. 

A 1-loop computation of all the graphs of the types in Eqs. (|A.3| )-( [A78D results in 
the following contributions to the effective parameters in Eqs. (|2.2j )- (|2.12| ) and to the 
wave functions in Eq. (ATI): 
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A. 2 Integrating out bosonic non-zero Matsubara modes 

The bosonic non-zero Matsubara modes are integrated out in the theory defined by the 
tree-level couplings in Eq. ( |A.2|) . Thus, there are loop contributions involving Higgses, 
right-handed stops and the gauge fields. 

The types of graphs and integrals appearing have been discussed in |33| (see also |J- 
JTTJ] ) . For completeness, let us briefly review the results here, as well. We use the 
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following notation: 
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where p& denotes the bosonic Matsubara momenta and a prime means that the zero 
Matsubara mode is omitted. Then the basic contributions are of the types 
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(A.33) 



Here the external momentum is q = (0, q). Note that in contrast to Sec. |A.1| , there 
are no subtractions to be made at 1-loop level: due to momentum conservation, all the 
internal lines have a non-zero Matsubara mode and are thus heavy. 

In most places, we keep only the dominant terms of the different qualitative types 
in the high temperature expansion in m 2 /(27rT) 2 , m 2 ,/ (27rT) 2 , m 2 2 / (27rT) 2 , m 2 / /(27rT) 2 . 
Note that to obtain these expressions, it is convenient to write the Higgs field propa- 
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gators as 



2 
TO 

, (HZHiC) = e mn e^j^ 2 . (A.34) 



We then obtain that 
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A. 3 Integrating out fermions 

The fermions to be integrated out, interacting with the Higgs and stop degrees of 
freedom, include the third generation quarks, the Higgsinos, and the SU(2) and SU(3) 
gauginos. The mass parameters of the latter are denoted by M 2 , M 3 . 
Concerning the integrals appearing, let us define 
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where pf denotes the fermionic Matsubara momenta. Then the basic graphs give, in 
the high temperature expansion, contributions of the form 



O-- 



£/[p 2 + M 2 ][(p + g) 2 + M 2 ] 
y + T ^(4M 2 2 + 4H 2 + 2g 2 )ln (4/i)2 



t4 



(A.54) 



1 



%[p 2 + Ml}[(p + q)i + 



16tt 2 



In 



~Pt 



Ml + 



3* 



(A.55) 



28 



i s [p 2 +Miw+w} (A - 56) 
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As the fermions appearing are Majorana particles, the Lorentz-structure of the con- 
tractions is 



(W a (p)Wp(q)) = 5(p-q) 



p 2 + M 2 ' 
-if 1 + M) 
p 2 + A 
[C-\ij> + M)U 



(W a ( P )W, {q) ) = S(p + g ) K=^^k , (A.57) 
<W a (p)W,(<z)> = *(p + ^ p2 + M2 , 

where the charge conjugation matrix C satisfies 

C = -C ] = -C T = -C~\ C^YC = -{Y) T , C~ l ^C = ( 7 5 ) T , (A.58) 

and W a (x) = Wp{x)Cpa. 

The final results for the couplings are (terms proportional to the gauge couplings 
come from loops involving gauginos, terms proportional to the Yukawa coupling from 
loops involving left-handed quarks): 



iZl = !_$!,„ MS!, SZ2 = ( 3 -f + 3k f)-LJ^l, (A .59) 



^ = ^(? + 4 « + l^l^ ln ¥)' (A.61) 
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^(T) = + 4(M| + l^^lnM!) + i^ 2 T 2 , (A.62) 



g 2 w *, (4/i) 2 



<5m 2 2 (T) = 3^M 2 /i* ln (A.63) 

/Xy / V 6 " ' 1071^ /iy 



^-3fM|„ m , i72 = ^( 2ft? + £,| + 3^).„(|!, (A.66) 
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For completeness, let us also write down the results for the CP-violating parameters 
in the limit that M 2 , are large compared with 7rT (this limit corresponds to the zero 
temperature contributions considered in p0[ ). This limit is of interest since it was seen 
in Eq. ( 5.17 ) that non- vanishing values of M 2 , \\x\ are preferred from the point of view of 
spontaneous CP- violation, and it is thus important to know when the high-temperature 
expansion used is reliable. We obtain 
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where <ip = d 4 ~ 2e p/ (2ir) 4: ~ 2e . Comparing with the expressions in Eqs. ( |A.59| )-( |A.71| ), we 
observe that the high-temperature expansion breaks down when {M 2 , fi} > IT. In the 
range 2T <{M 2 , < 4T, a numerical evaluation of the integrals appearing is needed, 
while at {M 2 , \fi\} > 4T, the values of the couplings saturate and the zero-temperature 
results in Eqs. ( |A.72| )-( [A.76D can be used. In Eqs. ( |A.72| )-( [A.76D we have shown also 
the direct integral representations which can easily be evaluated at any temperature 
(see, e.g., H). 



A. 4 Field redefinitions 



In Sees. |A.1HA.3| , the removal of various ultraviolet degrees of freedom from the theory 
resulted in contributions to the effective parameters characterizing the interactions of 
the infrared degrees of freedom. In addition, there were contributions to the field 
normalizations, 8Zi, SZ 2 , 5Zi 2 , bZy. It is convenient to make a redefinition of the fields 
so that these extra wave function terms disappear. This results in new contributions to 
the effective parameters, rendering them gauge and scale independent at 1-loop level. 

Summing the tree-level terms with those in Eq. ( A.l ), the kinetic terms are of the 
form 



C kia = (D?H i y(DfH 1 )(l + 5Z 1 ) + (DV'H 2 y(DV 1 H 2 )(l + 5Z. 



2j 



+ 



(D^H 1 )\D^H 2 )5Z 12 + H.c. 1 + {D!U)\DIU)(1 + 5Z V ). (A.77) 



One can now get rid of 8Z\, SZu by writing 



H x = (l-±5Z 1 )Ht w) -±5Z 12 Ht w \ (A.78) 

1 , T] -(new) A, 1 r y \ fV(now) 



H 2 = --5Z* 12 Hr w) + (l--SZ 2 )Hr w> , (A.79) 
U = (l--5Zu)u {ncw) . (A.80) 



2 

When the field redefinitions combine with the tree-level couplings in Eq. (|A.2|) , there 
are new contributions to the effective parameters as follows: 

5m\(T) = -m\8Z x - \m\ 2 (5Z 12 + 5Z* 12 ), (AM) 
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8m\ (T) = -m 2 2 5Z 2 - l -m\ 2 (5Z 12 + 5Z* 12 ) , ( A.82) 

5m 2 12 (T) = -^m 2 12 (5Z 1 + 5Z 2 ) - Urn 2 + m 2 2 )6Z 12 , (A.83) 

8m 2 v {T) = -mlrSZv, 5X V = ~g 2 s 5Zu, (A.84) 

5 7l = h 2 \fi\ 2 {5Z 1 + SZu) - h 2 (5Z* 12 A* t fi* + 8Z l2 A t jl), (A.85) 

5l2 = _ \A t \ 2 )(5Z 2 + 5^) - h 2 (5Z* 12 A* t fi* + 5Z 12 A t jl), (A.86) 
5 ll2 = -h 2 (l - \A t \ 2 - \fi\ 2 )5Z 12 - h 2 A* t il*(25Z u + 5Z X + 5Z 2 ), (A.87) 

S\ 1 = -^g 2 8Z 1 , 5X 2 = ~g 2 5Z 2 , (A.88) 

5A 3 = -\g\5Z, + SZ 2 ), 5\ 4 = \g\bZ x + 5Z 2 ). (A.89) 
There are no new contributions to A5...A7. 

A. 5 Gauge couplings and Co-integrations 

So far we have been considering couplings related to the Higgs fields and to the stop 
fields. Let us now review what happens with the gauge fields. 

We start by inspecting the spatial sector, where the plasma does not screen the gauge 
fields. What temperature does is that it changes the effective gauge coupling related 
to the static fields [ISl 



g% = g 2 0) ( 1 + -7^-3 /? g auge ( In =~ + ^) + ^scalar In =^ + Aermion In J~ ) 

= 9Tq (A? au g e /^scalar "I - /pfermion) If- T^T; (A. 90) 

where T is any reference temperature, and the standard contributions to (3 are, for 
the weak gauge coupling (N = 2), 

/? gauge = 22iV, (A.91) 

Scalar = -N H (Higgses) - 3iV g (squarks) - iV s (sleptons), (A. 92) 
Pfermion = -4N (gauginos) - 2iV H (Higgsinos) 

- 6 A" 9 (quarks) - 2 A" 9 (leptons) . (A.93) 

Here = 2 is the number of Higgs doublets and N g = 3 the number of generations. 
For the superpartners, often only one generation (squarks) or none at all (sleptons) 
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is taken to be effective. The expressions for the strong gauge coupling are completely 
analogous (see also 0). 

Consider then the temporal components of the gauge fields. They do get screened. 
Neglecting higher order (9(g 4 )-terms related to A ,Cq, the (tree- level) part of the La- 
grangian involving them is 

£a ,c = ipfiWiT^ 



+ ^g 2 A a A a (HlH 1 + HlH 2 )+glC^U a (T A T B UU;. (A.94) 
Here the covariant derivatives are in the adjoint representation and the Debye masses 



m 2 Ao = -g 2 T 2 , m 2 Co = -g 2 T 2 , (A.95) 



are (see, e.g., |48fl ) 

5 

2^ " 1 3" 

where we took into account gauginos and Higgsinos {Sm\ = (l/2)g 2 T 2 ), as well as 
gluinos (SmQ Q = (l/2)g|T 2 ), but only the 3rd generation ?7-squarks. 

Now, usually the fields A , C can be integrated out, since they have the mass scale 
~ gT, while the phase transition dynamics is related to the scale ~ g 2 T. However, we 
have argued in Sec. || that to study spontaneous CP-violation, one has to keep also 
the Higgs degrees of freedom with masses ~ gT in the action. Thus it is, in terms 
of the original 4d power counting, not strictly speaking consistent to integrate out 
some degrees of freedom with masses ~ gT and leave others in the action. However, 
we can introduce a different power counting within the 3d theory, in particular since 



numerically the Debye masses in Eq. ( |A.95| ) are relatively large. Moreover, they do 



not affect the CP-violating couplings at leading order. Thus it seems very reasonable 
to integrate out also A , C . The 1-loop results of that integration are |5|: 



247rm J 4 / V 167rmc . 

6m 2 (T) = 5m 2 2 (T) = -^g 2 Tm Ao , 5ml (T) = -^-g 2 s Tm Co , (A.97) 

l07T 37T 



SXu = (A.98) 
3d 87rm Co 

Q 4 r 7 - ' Q 4 r 7 - ' 

^ = 5^ = ---^ , SX 3 = --/ . (A.99) 



A. 6 Vacuum renormalization 

In the previous sections, we have derived an effective theory in terms of renormalized 
parameters in the MS scheme. In order to make connection to physics, the MS scheme 
parameters should, in turn, be expressed in terms of measurable quantities. For com- 
pleteness, we will discuss here some of the most important effects. A more detailed 
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discussion can be found, e.g., in [p| — [ITTfl . In this section, we take, for simplicity, the 
couplings A t , fi to be real: 

A; -> A t , li* -> Li. (A.100) 

Of course, these quantities could still be negative. For the case of complex couplings 
see, e.g. 
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Let us note first that the problem with the renormalization of the gauge couplings 
<7§, g 2 and the top quark Yukawa coupling h t is that the same parameters appear 
in many different vertices due to supersymmetry, while in the softly broken theory, 
this equivalence is lost beyond tree-level. Thus, it is not sufficient to compute three 
different physical observables to fix these three parameters appearing in different places. 
However, as the renormalization effects related to these parameters are not too essential 
for the problem discussed in this paper, we shall not discuss the fixing of g, h t , g$ in great 
detail: it is sufficient to remark that the effective finite temperature gauge couplings in 
the gauge sector are related to the zero temperature MS scheme couplings as discussed 
in Sec. |A.5|, and the top quark Yukawa coupling can be fixed at tree-level, 



,2^,2 



hfsm 2 (3 = (A.101) 



g m 

2 

Ml' 

where tan (3 = 1)2/1)1 is a parameter determined by the Higgs mass, see below. 

Let us, instead, discuss the values of the mass parameters mf(/2), m|(/2), mi 2 (fi), 
mfjlfl). For them, there are relatively large renormalization effects proportional to h 2 , 
and we would like to account for the dominant terms of such types, in particular those 
proportional to m 2 Q. 

The four running mass parameters can be fixed by computing the 1-loop pole masses 
m i> m A> m L m f °f the Z boson, the CP-odd Higgs particle, the lightest CP-even Higgs 
particle, and the lightest stop, respectively. For the dominant terms we are interested 
in, this computation can be simplified since the momentum dependence of the 2-point 
functions can be neglected: indeed, the momentum dependence results essentially in 
multiplicative terms proportional to p 2 ~ m 2 , where m 2 is the pole mass in ques- 
tion. These terms are suppressed with respect to the dominant additive ones, ~ rriq. 
Since the momentum dependence can be neglected, one can derive some of the 2-point 
functions from the effective potential. 

There is one further simplification. It turns out that the only large additive (~ rrin) 
contributions to the Z boson mass m 2 z come from the tadpole graphs which account 
for a shift in the location of the broken Higgs minimum (1)1,1)2)- In other words, in 
the one-particle-irreducible graphs, the terms proportional to ttiq cancel due to gauge 
invariance (for explicit expressions see, e.g., Eqs. (A.23)-(A.25) in 0). This means that 
if we denote by (v 1, v%) the location of the radiatively corrected broken minimum, then 
we can write, within the present approximation, m 2 z = ^ 2 (w 2 +w|)/4, where g 2 = g 2 +g' 2 . 
Then, all appearances of vi,v 2 can be expressed in terms of tan/3 = 1)2/1)1 and rriz, as 
gvi = 2mz cos/?, gv 2 = 2mz sin/5. 
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To solve for m\ (p,), m 2 (p,), m\ 2 {Ll), we now proceed as follows. First, we solve from 
the conditions dV(vi, v 2 )/dvi = 0, dV{y i, v 2 )/dv 2 = for m\ (/i), m 2 (p,) in terms of 
tan (3, m\, m\ 2 (ji). We feed these results into the mass matrix of the CP-odd Higgs 
(and of a Goldstone boson). Solving for the physical eigenvalue m\, we obtain m\ 2 (ji) 
function of tan /3, m 



Ai m Zi 



2 r. 



fx. With this expression, the results are known also for 
m\ m 2 (p,). It remains to eliminate tan/3: this is done by plugging the expressions 
found into the mass matrix of the CP-even Higgses, and thus finding the relation of 
tan [3 to the lightest CP-even Higgs mass m%. 

To proceed with this program, we write the Higgs doublets in component form as 



1 V2\-hl + ih\ 



1 ( h\ + ih\ 
1 ~ V2\h° 2 - ih\ 



(A.102) 



The tree-level potential ((h®) = V\, (h 2 ) = v 2 ) is then 

V tiee (v u v 2 ) = -m\v\ + ^m 2 2 v 2 2 + m\ 2 v x v 2 + ^f^l ~ «f) 2 - 

The h 2 -contributions of the 1-loop potential to dV/dvi can be read, 
Eqs. (A.29), (A.30) in M: 



(A.103) 
e.g., from 



^l-loop ( 



Vl,V 2 ) 



V 2 S h {fJ>), 



(A. 104) 



where 



S h (l2) 



h 2 

3- 1 



16tt 2 



h 2 
3- 1 



16tt 2 



11 + A t tan/3 
2 



2m 2 In 



top 



2 

m top 



— (m 



t L 



+ m L)( ln 



m^mi 



+ 1 



+A t (-A t + /i cot/3) (In 



In 



n 2 



mfmi 



+ 1 + 



m% + m 2 T mi 
T 1 In — - 



2 2 
m% — mj 



mn 



— (m^ — m t ~) In 



K 

mfmi 



+ 1 + 



In 



mn 



(A.105) 



(A.106) 



Here 



m 
rn 
m' n 



m Q + ™top 



+ c °s2/3 



is = m to P (A - //cot/3) 
o 1 



2 

2 1 
m,- = - 

* 2 



m? + m? + \/(m 2 



mi 



2 1 2 
77% + m 



top ) 



+ 4mf 



thR 



2 1 2 

m ( - + m,- — 



m~ — m~ ) 2 + 4m ~ 

tly l R/ l LR 



(A.107) 
(A.108) 
(A.109) 

(A.110) 



35 



Note that in the notation of |§, w s — > —h t fi,u s — > h t A t ,e s = d s — > 0, — > m~ , 
-8, - <• »*» - <»• ml + - m|, m»_ m?. In add.tion, no Higgsino and 
gaugino contributions were considered in |J, but this does not affect the leading h 2 
terms were are considering here (see, however, the discussion below). 
It then follows from the conditions d(V tree + V 1 ' 10 ^) / dvi = that 

^iCaO = —rnl 2 (p,) tan/3 — -m% cos2 f3 — Ch(ji), (A. Ill) 

m l(P') = —Tn^ifi) cot (3 + -m| cos2/3 — Shift)- (A. 112) 

The mass matrix for the CP-odd Higgs mass uia and for a neutral Goldstone boson 
is given by the fields h\ , h 2 , and is of the form 



M, 



( m\{p) + \~g\v 2 - vl) + UfM ml 2 (p) + Uf 2 (ft) \ 
dd 1, mUfi) + Uf 2 (fl) rnKfl) + \~g\v 2 - v 2 ) + Uj 2 {p) J ' 

(A.113) 

The expressions for H£j(p) can be obtained, e.g., from Eqs. (A.19)-(A.22) in || (the 
functions Fh vanish within the zero-momentum approximation considered, and the 
overall sign should be reversed): 

n£(/2) = C h {p) - taxtf A tt iD(fi), nf 2 (fc) = A tt iD{fl), (A.114) 
n4(/i) = S h (fl) - cotpAtfJiD(ji), (A.115) 



where 



D(p)=3 -M,(lnJ^ + l + 4±4ln^). (A.116) 
lo7T z V uifirii m;~ — m~ mf 



Plugging in the values of m 2 (fi),m 2 (fi) from Eqs. ( |A.111|) , ( |A.112| ), the CP-odd mass 
matrix becomes 

M f -[ml 2 (p) + A t fiD(p)}tan(3 m\ 2 {p) + A t /iD(p,) \ (AU7 s 

cp -° dd { m 2 12 (fi) + A tf iD(fi) -K 2 (/i) +A tf iD(fl)} cot/3 )' l ' ' } 

The (scale- independent) eigenvalues of Mcp_ dd are zero and m\, which finally gives the 
expression for mf 2 (fi), and after insertion into Eqs. ( |A.111| ), ( |A.112| ), for m\ (/2), m 2 (fx): 

m i2{P) = -~m 2 A sin2(3 -A t iiD(p,), (A. 118) 

ml(fi) = -m 2 A - -{m 2 A + m|) cos2/5 + \/i\ 2 D(fi), (A. 119) 



2^2 

ml(fl) = -m\ + X -{m\ + m%) cos2/? + \A t \ 2 D{ji) (A.120) 



+3 



h 



16tt 2 



2 \ tyy-t 2 



(m 2 + m 2 j) ( In — — h 1 ) + 2m? on In — + (mi — m?) In — - 

v mfiJii J 1 miirif 1 mf 
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There is the following useful observation to be made. If one expands Eqs. (| A . 1 1 8|) 
( |A.120| ) in terms of m 2 /mQ,m 2 op /mQ in the limit of large mg, then it turns out that 
the dominant terms containing mq cancel against the terms in Eqs. (|A.11|) -( |A.14|) . 
For example, summing the dominant 1-loop contributions to m 2 2 (T) from Eqs. ( |A.14[ ), 
( |A.39| ) together with Eq. ( |A.118|) , the scale dependence cancels and we obtain, up to 
terms of relative order m~ /wig, mf op /rriQ, 



Sm 2 



12 



(T) 



16tt 2 



A*/t* m 



In 4 + 1 

m\ 



4 m. 



top 



14 



rrir 



42 1 



.3. 

m 2 r, 



In 



mi 



(A.121) 



where we have denoted mf = m tory A t . The largest term here is that induced by the 

*LR " 

symmetry breaking, ~ A* //* (m top / nig) 2 . For m 2 (T),m 2 (T), in particular (where the 
result is similar), such terms are typically much smaller than the tree-level mass terms. 
Thus, after vacuum renormalization is taken into account, the terms proportional to 
h 2 /(lQiT 2 ) in m 2 (T),ml(T) can be neglected as a first approximation, and it is sufficient 
to consider the tree-level mass terms and the thermal T 2 -terms. 

By similar arguments, we can also easily derive the dominant contributions of the 
parameters M 2 , in the limit that M 2 , ^> my/- Indeed, in this limit the expressions 
for vacuum renormalization must cancel the large direct contributions in Eq. (|A.73 ), 
just as happened above for rriQ in the limit mq ^> m top . Thus, e.g., 



m 2 2 (/i) « m 2 2 (/i) 



(A.118) 



5m\ 2 {T) 



(A.122) 



Finally, let us connect tan/5 to the Higgs mass. The mass matrix for the CP-even 
Higgses m h , mu is given by the fields h^h®: 



Mi 



m 



CP-even 



2 Gu) 4 l~g 2 (3v 2 - v 2 ) + ng02) mUfi) - \~g 2 v lV2 4 U&Qi) 



m ? 2 (^) - 79 2 viv 2 4 nf 2 (/i) 



m 2 2 (fi) + l~g 2 (3v 2 - vf) + H*(fi) 



(A.123) 

where the nfj(/i)'s can be obtained from the second derivatives of the effective potential 
V(vi,V2), or directly from a diagrammatic computation in, e.g., Eqs. (A.12)-(A.18) 
in 0. Denoting 



5x 



12 



12 



16tt 2 
hi 



1 4 



m% 4 mf , mi 
T * In — - 



2 2 
m|, — m| 



TTlrj 



in 



top 



167r 2 m 2 ^, 



mi 



ln^. 

mi 



12- 



'16tt 5 



, mim^ 

ln^rr^. 



m 



top 



|/l t | 2 m 



top 



(m|, — m 2 ) 2 ' 



(A.124) 
(A.125) 
(A.126) 



37 



where again A t = A t — fi cot/3, we obtain 



nfiQx) = -H ^ - , n f 2 (^) = Atfi D(fi) -5, - A tf i5 2 , (A.127) 



-IA 

_Q L_ 

167T 2 



£>(£) - 5 



( m Q + m l) ( ln + + ("4 ~ m l) 



T mz) ln 



rrirp 



(A.128) 



Plugging these expressions together with those in Eqs. ( A.118| )-( |A.120| ) into 
Eq. ( |A.123| ), the explicit scale dependence cancels and the physical mass matrix re- 
maining is (see also fl5D| ) 



M, 



CP-even 

M n 

M12 

M 22 




M 12 \ 
/12 M 22 ) ' 
m| cos 2 /? -f-m^sin 2 /? + |/i| 2 5i, 
— (m| + m^) sin/? cos/3 

,2 • 2,3 i _2 „„„2 



2; 



m 2 z sin 2 /? + cos 2 /? + | A| 2 <5i + 2A t A t 5 2 + m 2 top S 3 . 



The lightest CP-even Higgs mass is then 



(A.129) 

(A.130) 
(A.131) 
(A.132) 



m. 



M u + M. 



22 



(M; 



11 



M 



22 J 



+ 4|M 12 | 2 



(A.133) 



The renormalization of mfj(jl) and its relation to the lightest stop mass (which, 
in the absence of mixing, is the right-handed stop mass mz ) can be handled with 
similar methods as the Higgs sector. However, we have already argued what the result 
will approximately be: there are large renormalization effects related to the integration 
out of Q, see Eq. ( A. 14 ), but when vacuum renormalization is taken into account, 
then these effectively cancel. Thus, it is enough for our purposes to ignore the effects 
in Eq. ( |A.14[) and replace mfj with the tree-level value, which is obtained from the 
expression for the lightest stop mass: 



4 



Hop 



1 - (A t - fi cot (3f 



(A.134) 



A. 7 Summary: the effective theory 

In this section we collect together the results of the previous sections. For the quartic 
couplings, we will work at the accuracy of the vacuum renormalization discussed in 
Sec. |A.6| . Thus, we ignore all the contributions of order 0(g 4 ) to the quartic Higgs 



couplings, contributions of order 0(g$) to the quartic stop coupling, and contributions 
of order 0(hf) to the top Yukawa coupling squared appearing in the quartic scalar 
interactions between stops and Higgses. In other words, we keep only contributions for 
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which the scale dependence is cancelled within the accuracy of vacuum renormalization 
discussed in Sec. |A.6| . It should be noted, though, that for the quartic stop coupling, 
neglecting effects of order 0(gg) and keeping effects of order O(hf) is numerically not 
a very useful approximation scheme. 

We also make an expansion in m./mL m 2 /(nT) 2 , where m 2 denotes the mass pa- 
rameters other than mg. In general, we keep the dominant effects of the different 
qualitative types appearing. Note that \A t \ can, in principle, be larger than since 
it is not subject to constraints from the high-temperature expansion. The Debye masses 
appearing are given in Eq. ( |A.95|) . 

For the mass parameters m^iT),™ 2 ,^),™^^), we only keep the leading terms 
below, which should already give a reasonable approximation as explained after 
Eq. (|A.121|) . More precise results can be obtained by using the full expressions in 
Eqs. ( POT9l ),( P02q ) (see also @H[lTJ). 

For m 2 2 (T), which affects directly spontaneous CP- violation, we display a more 
precise expression, with the dominant additive corrections of various qualitative types 
included. Multiplicative corrections to m 2 A sin2/3 have been neglected. It should also 
be noted that the high temperature expansion for the gaugino-Higgsino contribution is 
reliable only when the expression in the square brackets on the first row in Eq. ( |A.139|) 
is negative: without a high temperature expansion, this term is 



2 ,_. 2 M 2 /i 



where 



(A.135) 



I f ( y ) = / dxx 2 = = — . (A.136) 

If M 2 , \fi\ are large (> ttT), then the term in Eq. ( A.135 ) vanishes, but also the contribu- 



tions of these degrees of freedom to the T 2 -terms in the mass parameters, Eqs. ( |A.61| )- 
( |A.64j ), and to the Debye masses, Eq. ( |A.95|) , should be left out. 
Note that the square bracket on the second row in Eq. ( |A.139| ) is also negative, so that 



the whole term is positive for sign (A*/}*) > 0; see Eq. ( |A.121|) for an approximation. 



Finally, let us point out that we write here the couplings in 4d units, so that there 
is an overall factor Jq dr = l/T in the action. This factor can, as usual, be defined 
away by a rescaling of the fields, and the result is that the quartic couplings of the 3d 
action get simply multiplied by T. 

With these conventions and approximations, we obtain that 8 

1 X 

m l( T ) = 7^4 - -(m 2 A + m 2 z )cos2p 



3 Note that as can be seen from the hf-tevm in Eq. ( A.146| ), there is an error in Eq. (3.7) of the 



latter of [jl2|: there should be no extra factor of e 3 / 4 inside the logarithm. (This does not affect any 
of the conclusions in O] concerning the non-perturbative effects.) 
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(A.140) 



^(41n^ + 3K| 2 + 3|/i| : 
+ Qhl A n ^ + 2 \A t \ 2 In ^ + |/i| 2 - (|i t | 2 + |/i| 2 ) 2 In , (A.141) 

,2 



(A.142) 
(A.143) 

(A.144) 
(A.145) 

(A.146) 
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In these formulas, all the scale dependence has cancelled. 
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